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The vortex dynamics in a low temperature superconductor deposited on top of a rectangular 
array of micrometer size permalloy triangles is investigated experimentally. The rectangular unit 
cell is such that neighboring triangles physically touch each other along one direction. This design 
stabilizes remanent states which differ from the magnetic vortex state typical of individual non- 
interacting triangles. Magnetic Force Microscopy images have revealed that the magnetic landscape 
of the template can be switched to an ordered configuration after magnetizing the sample with an 
in-plane field. The ordered phase exhibits a broad flux flow regime with relatively low critical current 
and a highly anisotropic response. This behavior is caused by the spontaneous formation of two 
separated rows of vortices and antivortices along each line of connected triangles. The existence of a 
clear flux flow regime even for zero external field supports this interpretation. The density of induced 
vortex-antivortex pairs is directly obtained using a high frequency measurement technique which 
allows us to resolve the discrete motion of vortices. Strikingly, the presence of vortex-antivortex 
rows gives rise to a self organized synchronized motion of vortices which manifests itself as field 
independent Shapiro steps in the current- voltage characteristics. 

PACS numbers: 74.78.Na 74.78.Fk 74.25.Dw 74.25.0p 



I. INTRODUCTION 

The dynamics and pinning of vortices in typc-II su- 
perconductors is arguably one of the most widely inves- 
tigated phenomena in the field of superconductivity. Re- 
cently, particular attention has been devoted to the pos- 
sibility of controlling the strength of the vortex pinning 
using different magnetic templatesii^. Within the Lon- 
don formalism, the interaction between a vortex line and 
a finite size permanent magnet is given by^. 



«7p(R) 



m(r) •B„(R-r) 



(1) 



dot 



where the integration is carried out over the volume of 
the ferromagnet, R indicates the position of the vortex 
line, m(r) its spatial dependent magnetic moment, and 
B„ is the field generated by the vortex line. This equa- 
tion shows that the pinning potential not only depends 
on the size of the dots^i^ but also on their exact magnetic 
stated. In particular, if the superconducting penetration 
depth A is on the order of the lateral size of the dot, a 
weak average pinning is expected when the magnetic dot 
contains multiple domains. In contrast to that, a maxi- 
mum pinning should be achieved for single domain struc- 
tures. In order to retain these properties in big arrays of 
magnetic particles it is necessary to ensure a minimum 
distance between them in such a way that the field pro- 
files of neighboring elements do not overlap and can be 



resolved at scales of the magnetic penetration length A. 

Notice that according to Eq. ([T|) it is possible to increase 
strongly the pinning by just growing continuously the 
magnetic moment of the dots. This growth is limited by 
the fact that above certain magnetic moment, a vortex- 
antivortex pair will be induced by the magnetic element 
with the consequent change in the pinning properties^. 
For instance, dots with perpendicular magnetic moment 
induce a vortex on top of the dots whereas the accom- 
panying antivortex is located at interstitial positions, in 
between the dots. As a result, vortices and antivortices 
feel a different pinning potential which leads to a dis- 
tinct dynamic response^.. Since the weakly pinned in- 
terstitial antivortices limit the maximum critical current 
jmax ^ ^jjgjj. annihilation with an external field produces a 
shift of the field position of the J™°^ towards H = nHi , 
where n is the number of vortex-antivortex pairs gener- 
ated by a single magnetic element and Hi is the field 
corresponding to one external vortex per unit cell. Since 
Hi is a geometrically predefined parameter, the field loci 
of the J™*^^ is a direct measurement of the number n of 
bounded pairs. Unfortunately, in the case of micromag- 
ncts with in-plane magnetic moment, the field polarity 
symmetry of the system conceals this information from 
dc electro-transport experiments. 

In this work, we introduce alternative methods to ac- 
curately determine the number n of vortex-antivortex 
pairs induced by the magnetic template, based on the 
ac-dynamic response of the hybrid system. Indeed, the 
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fact that the critical force to bring the vortex and an- 
tivortex closer differs from that necessary to pull them 
aparli^, gives rise to a clear rectification effect when the 
system is submitted to a zero average ac-signal. This rec- 
tified voltage is maximum at zero external field^ and pro- 
gressively decreases with increasing field until all vortices 
(or antivortices) are compensated. At RF frequencies the 
mode locking or synchronized motion of vortices produces 
distinct voltage steps in the current-voltage characteris- 
tics with separation determined by the number of moving 
vortices. Both techniques indicate that for the particular 
system here investigated, the magnetic template gener- 
ates one single vortex-antivortex pair per unit cell. 

The crucial conditions for the success of these meth- 
ods are (i) a periodic landscape, i.e. all vortices and 
antivortices feeling the same environment, and (ii) a free 
flux flow motion of vortices. Although the predefined 
topological landscape is by definition periodic, the mag- 
netic landscape for the as-grown state consists of a highly 
disorder distribution of magnetic poles. This disordered 
state can be switched to a perfectly ordered magnetic 
lattice by magnetizing the sample with an in-plane field. 
The order-disorder transition, already reported for other 
S/F systems^, produces a profound modification in the 
vortex dynamics. More specifically, in the disordered 
state, non-linear current- voltage characteristics with high 
critical currents are observed, whereas a very broad flux 
flow regime accompanied by a sharp decrease of the crit- 
ical currents is obtained in the ordered state. This be- 
havior can be attributed to the reduction of the average 
hopping distance determined by the separation between 
neighboring magnetic poles of identical polarity. In ad- 
dition, the presence of a flux flow regime at zero applied 
field reinforces the idea that the investigated magnetic 
structure is able to induce vortex-antivortex (v-av) pairs. 



II. SAMPLE PREPARATION AND 
CHARACTERIZATION 

The sample investigated consists of a periodic array of 
permalloy (NigoFc2o) equilateral triangles of 1 fiui side 
and a thickness of 25 nm. Along one of the principal di- 
rections of the array (x-axis) the separation between the 
triangles is 1 fim whereas in the perpendicular direction 
(y-axis) the triangles are touching [see Fig. [TJa)] . This 
geometry gives rise to two well defined directions of the 
structure, parallel (y) and perpendicular (x) to the hue 
connecting nearby triangles, as indicated in Fig. [IJa). 
A 50 nm thick superconducting Al bridge deposited on 
top of the magnetic template also covers the unpattcrned 
substrate [see Fig. (Hb)] thus allowing us to compare di- 
rectly patterned vs reference plain film and identify the 
influence of the magnetic triangles on the superconduc- 
tor. Both, the array of triangles and the transport bridge 
are aligned at submicrometer scale and were prepared by 
standard lithographic techniques. The Al plain film has 
a critical temperature at zero field Tco = 1.315 K, with 



residual resistivity ratio i?(300K)/i?(4.2K) ^ 2, and su- 
perconducting coherence length ^(0)w 128 nm. 

It is a crucial point that neighboring triangles over- 
lap in one direction in order to eventually achieve two 
well distinguished phases: disordered in the as-grown 
state and ordered in the magnetized state. Indeed, Mag- 
netic Force Microscopy (MFM) images shown in Fig. [TJc) 
for non touching triangles demonstrate that the rema- 
nent state, irrespective of the magnetic history of the 
triangles, is always a vortex state, in agreement with 
previous reports on individual micrometer size magnetic 
structuresiiii^. In contrast to that, if there is a phys- 
ical contact between neighboring triangles a new set of 
possible states can be found. Fig. [T](d) shows a MFM 
image of the sample investigated in the as-grown state. 
Here a clearly disordered arrangement of the, so called, 
buckle statesH is observed. Interestingly, after magne- 
tizing the sample with an homogeneous in-plane field 
perpendicular to the chain of triangles [Fig. [IJe)] , these 
buckle states can be set in a perfectly ordered arrange- 
ment. The in-plane distribution of magnetic moments 
calculated by OOMMF micromagnetic simulations.ii is 
indicated in Fig. [ijf) with black arrows whereas the z- 
component of the stray field (Bz) is shown in Fig. [Ijg). 
In all cases the magnetic state of each triangle is charac- 
terized by a strong out of plane magnetic pole which is 
located at one of the two bottom corners of the Py trian- 
gles and a weaker magnetic pole, with opposite polarity, 
located at the top corner of the triangle [Fig. [TJe)] . It 
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FIG. 1: (color online) (a) schematic drawing of the array of 
Py triangles. The direction of the external field, the current, 
and the relative orientation of the array is clearly indicated 
with arrows, (b) Schematic representation of the Al trans- 
port bridge with Hall voltage contacts on the patterned re- 
gion. Magnetic Force Microscopy pictures obtained at rema- 
nence of physically separated magnetic triangles (c) , touching 
magnetic triangles in the virgin state (d), and the triangles 
magnetized perpendicular to the chain of triangles (e). Panels 
(f) and (g) show the local magnetic moment distribution and 
the z-component of the stray field, respectively, calculated by 
OOMMF simulationsi^. 
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is precisely this magnetic pole which can be controlled 
to lie on either side of the triangle by simply changing 
the direction of magnetization. Notice that this magnetic 
monopole-like structure breaks the field polarity symme- 
try of the S/F hybrid system. The question now arises 
as whether the distribution of the magnetic dipoles (or- 
dered versus disordered) has any effect on the properties 
of the adjacent superconducting film. 



III. NUCLEATION OF SUPERCONDUCTIVITY 



Let us first analyze the nucleation of the supercon- 
ducting order parameter for the two different magnetic 
states of the hybrid system. Fig. [5] summarizes the ob- 
tained superconductor/normal metal phase boundaries 
for a dc-currcnt of 100 /iA and using 90% criterion of the 
normal state resistance i?„ for different magnetic states 
and current orientations. Higher criteria exhibit simi- 
lar results. For the as-grown state (open circles) a pro- 
nounced decrease of the critical temperature at low fields 
with respect to the plain film is observed. A further 
suppression of Tc is obtained after the sample has been 
magnetized with an in-plane field parallel to the x direc- 
tion (solid triangles). The fact that the phase boundary 
of the as-grown state lies in between those corresponding 
to the plain film and the fully polarized buckle state, sug- 
gests that in the virgin sample, randomly oriented buckle 
states coexist with magnetic vortex states. This has been 
indeed confirmed by the MFM images taken at different 
locations of the pattern. 

It is worth noticing that for the magnetized sample the 
Tc{H) boundary exhibits a parabolic background similar 
to that reported in earlier invcstigation o^^i^^i^^ . This be- 
havior can be ascribed to the change of dimensionality 
when ^(r) exceeds the size w of the regions where su- 
perconductivity first nucleates. Within this regime the 
phase boundary can be approximated by Tc{H)/Tc{0) = 
1 — {aH)"^, with a = ^(0)7rw/>/3(^o- By using this ex- 
pression to fit the data in Fig. [5] (solid line) we found 
w ^ lA ^m which is consistent with the largest separa- 
tion between neighboring magnetic poles. This behavior 
suggests that for the magnetized sample the stray fields 
confine the superconducting condensate in between the 
magnetic poles and consequently the superconducting or- 
der parameter will be modulated into a sequence of quasi 
one-dimensional stripes. This prediction is in agreement 
with complementary transport measurements for two or- 
thogonal orientations of the applied current, perpendic- 
ular to the chain of triangles (solid triangular symbols) 
and parallel to them (open triangular symbols) as shown 
in Fig. [2] indicating that the onset of superconductivity 
is independent on the direction of the applied current. 
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FIG. 2: Superconducting/normal Pliase boundary as deter- 
mined by 90 % resistance criterion of the normal resistance 
Rn and dc-current of 100 for the as-grown state (open cir- 
cles) and different magnetized states and current orientations. 
The dashed line shows the corresponding transition lines for 
the co-evaporated Al plain film. The solid line is a fit assum- 
ing a lateral confinement of the superconducting condensate 
(see text). 



IV. VORTEX DYNAMICS IN THE ORDERED 
AND DISORDERED PHASES 

In the previous section we showed that by switching 
from the as-grown (disordered phase) to the magnetized 
(ordered phase) state, a proliferation of magnetic poles 
takes place which in turn gives rise to a decrease of Tc- In 
addition, the dimensional crossover from a 2D disordered 
phase to a ID ordered phase has little or no influence 
on the nucleation temperature. However, deeper into 
the superconducting state where the electric response is 
dominated by the vortex dynamics, the one dimensional 
nature of the field modulation should become more ap- 
parent. 

In order to address this issue we have measured the 
voltage-current characteristics [T^(/)] with the bias cur- 
rent applied in the a;-direction at = ■^Hi and T = 1.1 
K for the two different magnetic states: as-grown (open 
circles) and magnetized (solid squares) states as shown 
in Fig. [3] (a). The most obvious feature in this Figure 
is the enormous difference (more than an order of mag- 
nitude) between the critical current of the magnetized 
state (24 /iA) with respect to the demagnetized state (490 
/iA). Furthermore, the vortex dynamics is severely mod- 
ified after magnetizing the sample. Indeed, whereas in 
the as-grown state there is a narrow window of currents 
exhibiting a highly non-linear V{I) dependence, where 
the vortex lattice transit from pinned phase to the nor- 
mal state, in the magnetized sample a broad fiux-fiow 
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regime is observed immediately above the critical cur- 
rent. Fig. [2Ib) and (c) summarize the field dependence 
of the different dynamical phases for the as grown sample 
[panel (b)] and after magnetization [panel (c)]. The onset 
of vortex motion is delimited by a resistance criterion of 
10% R„ while a 90% R„ criterion is used to indicate the 
transition to the normal state. 



fectly aligned and strongly overlap, thus creating an easy 
channel for flux motion. In contrast to that, in the disor- 
dered phase, the average distance between two neighbor- 
ing magnetic pinning centers of equal polarity is larger 
thus impeding the easy motion of vortices. From now on 
we will focus on the dynamic response of the system in 
the ordered state. 
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V. EVIDENCE OF INDUCED V-AV PAIRS 

A rough criterion to determine whether the stray field 
emanating from the magnetic structure is strong enough 
to induce v-av pairs can be obtained by comparing it to 
the penetration field Hp resulting from Bean-Livingston 
barriersi^. The highest barrier possible is limited by the 
thermodynamic critical field Hp = '^q/ {2\/2-n\£^). Using 
^ = 128 nm and A = 150 nm we obtain Hp ^ 12 mT 
which lies well below the field intensity 100 mT) pro- 
duced by the magnetic template at the superconducting 
film. Based on this estimation we can safely conclude 
that the magnetic structure is able to induce, at least, 
one vortex-antivortex per unit cell. 




FIG. 4: (color online) Schematic drawing of the induced vor- 
tex lattice in the case of no external applied field (a) and an 
applied field equal to \Hi (b). Red (green) circles represent 
positive (negative) vortices. 



FIG. 3: (color online) (a) Voltage-current V{I) characteristic 
with current applied in the 2;-direction, at H — ^Hi and a 
temperature of 1.1 K for the as-grown and the magnetized 
states. The lower panels show the magnetic field dependence 
of the critical current and the transition to the normal state as 
determined by criterion of 10% R„ and 90% Rn, respectively, 
for the as-grown (b) and magnetized (c) states. 

Notice that the presence of matching features, which 
require two-dimensional phase coherence, together with 
the fact that in the magnetized state a finite critical cur- 
rent is needed to dcpin the vortices, rule out the pos- 
sibility of having normal metal stripes in the hybrid 
sample. Instead, the observed difference between both 
phases can be attributed to the formation of easily mov- 
ing channels of vortices and antivortices created by the 
magnetic structure. Indeed, in the ordered phase the 
potential wells created by the magnetic dipolcs are pcr- 



This finding suggests that the ordered phase consists 
of two well defined channels with opposite vortex polarity 
as schematically shown in Fig. |4l[a). Since the magnetic 
template breaks the field polarity symmetry, in princi- 
ple vortices and antivortices should experience a differ- 
ent depinning current. Moreover, as shown in Fig. [3I[b), 
under these circumstances an external field will compen- 
sate vortices in the channel with opposite polarity (thus 
introducing disorder in this particular channel and in- 
creasing the average hopping distance) and as a result 
the high vortex mobility will be suppressed. In contrast 
to that, the row with vortices of opposite polarity would 
remain unaltered. This scenario predicts two anomalous 
behaviors: the existence of a flux-flow response at zero 
field and a field independent flux flow regime as long as 
none of the two rows of vortices is fully compensated 
by the external field. Both predictions are in agreement 
with the field evolution of the V{I) measurements done 
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FIG. 5; (a) The field dependence of the flux flow resistance 
for the magnetized sample at T = 1.1 K. The straight lines 
indicate the expected dependence according to free flux flow 
model, (b) Critical current as a function of magnetic field at 
a temperature of 1.1 K for both current directions. A clear 
anisotropic response is observed. 



at 1.1 K. The slope of the linear regime, proportional to 
the amount of particles contributing to the easy flow, is 
plotted as a function of field in Fig. ^a). In compar- 
ison to the slope expected for the free flux flow model 
Rff = RnH/Hc2 [straight line in Fig. [5l^a)], where i?„ is 
the normal state resistance and Hc2 is the upper critical 
field, a weak field dependence of the flux flow slope is 
seen at low field values. 

Further evidence in favor of the existence of channels 
populated with v-av pairs comes from the measurements 
of the critical current for currents applied parallel (Jj,) 
and perpendicular {Jx) to the channels of v-av pairs. 
These results are shown in Fig.[S{b), with a critical cur- 
rent determined using a criterion of 100 /iV. This figure 
shows a strong in-plane anisotropy 7 = Jx/Jy induced 
by the magnetic template which results from the easier 
motion of vortices along the line of magnetic poles than 
perpendicular to them. This anisotropy amounts to a 
factor of 5 at non matching fields and is minimum at the 
commensurability fields. 



VI. DETERMINATION OF THE V-AV PAIRS 
DENSITY AT LOW FREQUENCIES: RATCHETS 

The scenario described in the previous section where 
rows of vortices and antivortices generated by the mag- 
netic template coexist, makes our system similar to the 
array of elongated magnetic bars studied by de Souza 
Silva et al^. This system has been modeled within the 
London approximation by Carneiro^ for the case where 
no v-av pairs are induced, and more recently by Lima 
and de Souza Silvaii^ assuming that v-av pairs emanate 
from the dipoles. Irrespective of whether the magnetic 



template is able to induce v-av pairs or not, both situ- 
ations lead to a rectification of the net vortex motion v 
when a zero average ac-force is applied parallel to the 
dipolar moment (i.e. along the line connecting a vor- 
tex with its corresponding antivortex)2£. However, clear 
differences in the field dependence of the average veloc- 
ity are expected. Indeed, if v-av pairs are present then 
(i) a maximum rectification is expected at zero field, (ii) 
unlike standard ratchet signal, the sign of the average 
voltage readout should be field polarity independent, (iii) 
as vortices (antivortices) are progressively compensated 
by increasing the external negative (positive) field, the 
rectified signal should decrease. 




FIG. 6: (color online) Rectified voltage Vdc in the 
temperature-field (T — H) plane for the magnetized sample 
subjected to a 1 kHz sinusoidal ac current in the y direction 
of amplitude 100 ^A. 



The measured rectification effect, shown in FiglHl con- 
firms all the above mentioned facts and indicates that 
the v-av pairs are more easily brought together than sepa- 
rated. Furthermore, if the applied field is increased above 
the first matching field a sign reversal occurs. This sign 
reversal at Hi can not be explained assuming more than 
one v-av pair created by the magnetic template since in 
this case the v-av annihilation process would be unlikely 
interrupted by the partial compensation of one row of 
vortices. In contrast to that, if there is a single v-av pair 
per unit cell, once a whole line is compensated by the ex- 
ternal field, a dramatic change in the rectification proper- 
ties is expected. Although this finding seems to indicate 
that the magnetic structure is able to induce a single v- 
av per triangular element, a more convincing evidence of 
this prediction can be obtained only by monitoring the 
individual motion of vortices at higher frequencies. 
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VII. DETERMINATION OF THE V-AV PAIR 
DENSITY AT HIGH FREQUENCIES: MODE 
LOCKING 

The motion of vortices in a periodic potential gives 
rise to a variety of interesting dynamical effects such 
as the generation of rf radiationSi and synchroniza- 
tion/resonance effects of the vortex Iattice22ii^i2i when 
a rf current 7^/ is superimposed onto a dc drive Idc- 
The latter effect gives rise to a series of plateaus in the 
V{I) measurements very similar to the Shapiro steps ob- 
served in Joshepson junctions^^. The relevant timescale 
of these processes is determined by the washboard fre- 
quency, f = v/d with d the period of the potential and 
V the average velocity of the particles. The investigation 
of the above mentioned phenomena can be used to ob- 
tain additional information on the dynamics of the vortex 
lattice^, the periodicity of the underlying potential and 
the density of particles in the system^!. In other words, 
high frequency experiments provide further insight into 
the microscopic motion of vortices not available with the 
averaged response presented in previous sections. 

The synchronous motion of vortices resulting from the 
resonant condition / = v/d^ known as mode locking ef- 
fect, allows one to obtain the exact number of v-av pairs 
created by the magnetic template. Indeed, the value of 
the average velocity is given by: 



V ~ cfd 



(2) 



with d the periodicity of the lattice in the direction of 
motion, / the frequency of the rf current and c the num- 
ber of cells that a vortex has traveled during a rf cycle. 
As we mentioned above, the fingerprint of this type of 
resonant motion is the presence of clear plateaus in the 
I{V) characteristics separated by voltage steps, 



Vc = Nk 



h V 
2ed 



Nk-^cf, 
2e 



(3) 



with k the number of vortices in each unit cell and N 
the number of rows between two voltage contacts. Mea- 
suring the voltage separation between consequtive steps 
and using Eq. ([3]) we can estimate the density of vortices 
moving coherently in the flux flow regime. 

Fig.lDJa) shows a V{I) scan for H =0.5 mT and T =1.1 
K with a small ac excitation Irf << Idc of 100 MHz su- 
perimposed. The presence of Shapiro steps in the linear 
part of the V{I) curve becomes more apparent by plot- 
ting V versus dl/dV [see Fig. [Tljb)]. These steps repre- 
sent a clear indication of the coherent flow of vortices. 

In Fig. mb) we assign the present peaks to values of 
Vc with c ranging from 1 to 4 (Vi = 63fiV, V2 = 127 fiV, 
V3 = IdA^V and V4 = 261nV). From their linear de- 
pendence on c, shown in the inset of Fig. [8l we calculate 
a value of fc « 1.06. An independent estimation can be 
obtained by tracking the frequency dependence of each 
voltage plateau Vc- This is shown in Fig. [5] for the first 



500 
400 
300 

> 

^ 200 



100 - 



(a) ■ ■ y 

/ 
/ 




' (b) 
- 


/ 






" / " 

















0.0 



0.2 

/ (mA) 



0.4 



V. 



^2 

V, 



1 2 

dl/dV 



FIG. 7: (a) V(I) measurement (filled symbols) at = 0.5 mT 
and T = 1.1 K with a small ac excitation of 100 MHz. Panel 
(b) shows a dl /dV , the discrete peaks are labeled according 
to Eq.®. 



two steps. The linear fit using Eq.Q gives fc w 1.1 and 
1.18 for Vi and V2, respectively. 

It is important to notice that, in general, two in- 
gredients are necessary for the presence of mode lock- 
ing effects. Firstly, the existence of a flux flow regime, 
which ensures a narrow distribution of average vortex 
velocities'^. Secondly, there should be a periodic pertur- 
bation of the pinning potential, otherwise the distortion 
of the vortex lattice becomes too strong and its coherent 
motion is destroyed. In principle, this second condition 
is only fulfilled at the matching field. Hi — 1.195 mT. 
This seems to be in contradiction with the fact that clear 
Shapiro steps are present at 0.5 mT, a field which would 
corresponds to 0.42 vortices per unit cell. This apparent 
paradox can be solved invoking once again the presence 
of v-av pairs induced by the triangles. In this case, the 
applied field will annihilate the vortices in the channel 
with opposite polarity, leading to an increase of their 
critical current and consequently suppressing their con- 
tribution to the coherent motion. On the other hand, the 
row with equal polarity will be unaltered by the applied 
magnetic field and therefore should exhibit a region of 
field independent Shapiro steps. In other words, the in- 
coming vortices generated by the external magnetic field 
play no role in the coherent oscillatory motion of vortices 
thus giving rise to a constant voltage separation between 
two consecutive Shapiro steps. 

In order to check this hypotheses of self organized com- 
mensurability we measured the field dependence of the 
Shapiro steps, which are better visualized as a peak in 
the dl/dV as previously shown in Fig. [Tfb). The results 
of these measurements for T =1.1 K are summarized in 
Fig. [5] as a contour plot in the field- voltage plane, with 
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FIG. 8: Frequency dependence of Vi (circular symbols) and 
V2 (square symbols) at a field of 0.5 mT and temperature of 
1.1 K. Inset: Vi -> V4 measured at 100 MHz. 




FIG. 9: (color online) Magnetic field dependence of dl/dV 
versus V at a temperature of 1.1 K and an applied rf frequency 
of 100 MHz. 



maximum intensity representing the exact location of the 
Vc peaks. For magnetic field values ranging from —Hi to 
+Hi a clear field independence peak position is observed. 
This result confirms that the amount of particles, moving 
coherently, is independent of the applied external field. 



VIII. DISCRETE MOTION OF THE 
VORTEX- ANTIVORTEX LATTICE 

It has been predicted by molecular dynamics simu- 
lations that a vortex lattice in a periodic landscape of 
asymmetric pinning sites can unveil the discrete hopping 
of vortices when the system is excited with a zero mean 
rf drive^Si^i^. So far, most of the experimental reports 
evidence a vortex rectification effect obtained in the low 
frequencies limit which results from an averaging in time 
and are not able to resolve the individual motion^^iS. In 
this section, we will try to bridge this gap by constructing 
an artificial vortex rectifier resulting from a dc tilt, which 
breaks the inversion symmetry of our system^!. Addi- 
tionally, this tilt allows us to probe the discrete motion 



of the vortex lattice without the need of a time resolved 
signal read out. 
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FIG. 10: DC voltage as a function of ac amplitude (/ = 100 
MHz) with a dc offset of 10 applied. The indexing (a, b) 
indicates the number of steps forward (a) and backward (b). 

Fig. [10] shows the dc voltage as a function of the 
applied ac drive with an additional dc current on top 
{Idc = 10 /iA). Due to the applied dc tilt the critical 
current, needed to depin vortices, is smaller for positive 
than for negative currents (/pos < Ineg)- Depending on 
the strength of the applied ac amplitude lac we can de- 
termine three different dynamic regimes. For small ac 
excitations {lac < I pos < Ineg) the driving force is un- 
able to depin the vortices and as a result the dc voltage is 
zero. By increasing the ac amplitude above Ipos the vor- 
tex lattice is set in the flux flow regime and mode locking 
effects result in an oscillatory dependence of the dc volt- 
age as a function of the ac amplitude. Finally at high ac 
amplitudes the normal state is reached and Vdc = Rnldc- 
Each maxima and minima in dc voltage can be attributed 
to the discrete motion of the vortex lattice in the positive 
and negative cycle of the driving amplitude. To identify 
the motion, we introduce the indexing (a, b) in Fig. 1101 
with a and b the number of steps in the positive and neg- 
ative direction respectively. At lac just above Ipos the 
driving force is strong enough to move the vortex lattice 
one unit cell in the positive direction while no motion 
occurs in the negative cycle of the ac drive and a max- 
ima in the dc voltage output occurs. If the ac amplitude 
passes Ineg a single back and forth step of the vortex 
lattice occurs and as such a strong reduction of the dc 
voltage is measured. This difference between back and 
forth motion of the vortex lattice, due to a small dc tilt, 
results in a discrete stepping procedure up to the sixth 
order. 

Fig. [11] shows the dependence of the dc voltage, nor- 
malized to Vi, as a function of the ac drive and the dc 
tilt. At small dc tilts (below w 100 /iA) the behavior 
is exactly the same as shown in Fig. 1101 In this region, 
the applied dc current creates a small asymmetry and 
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FIG. 11: (color online) lac and Idc dependence ol the dc volt- 
age at a magnetic field value of 1.195 mT and a rf modulation 
of 150 MHz 

the vortex lattice, depending on the exact value of the 
ac drive, can either move one step extra in the positive 
direction (maxima in dc voltage) or an equal amount of 
steps in both directions (minima in dc voltage). Aroimd 
100 /iA the oscihations vanish (white region) and the dc 
voltage readout amounts Vi . The reason for this behav- 
ior is a subtle balance between the ac drive and the dc 
tilt. Exactly at the point were the ac amplitude is strong 
enough to move the vortex lattice in the 'hard' direction 
(i.e. against the dc tilt) it is also sufficiently strong to 
move the vortex lattice an additional step in the 'easy' 
direction. Therefore the overall motion of the vortex lat- 
tice will be independent on the applied ac drive and equal 
to 1. This subtle balance is destroyed in favor of 'easy' 
flow if the applied dc current is increased above 100 /iA. 
Again the dc voltage readout (1 < Vdc/Vi < 2) indicates 
that the vortex lattice moves two steps (maxima in dc 
voltage) or one step (minima in dc voltage) extra in the 
positive direction. By increasing the dc tilt even further 
we could even see a difference of 4 unit cells between the 
motion in the positive and negative direction. 

IX. CONCLUSION 

We have demonstrated that the particular geometry 
of the two-fold symmetric array of magnetic triangles 



can induce an order-disorder transition as the system 
switches from the as-grown state to the magnetized state. 
The disordered phase is characterized by a random distri- 
bution of two opposite chirality magnetic buckle states, 
whereas in the order state all triangles are in a well de- 
fined unique magnetic buckle state. The latter phase ex- 
hibits unidimensional lines of strong stray fields which 
represent channels for easy flux motion resulting in a 
clear anisotropy of the critical current. A quantitative 
estimation shows that the current needed to induce inter- 
channels hopping is about 5 times larger than that needed 
for setting intra-channel motion. It is important to stress 
that even though we have presented clear evidence of the 
influence of the ordering of magnetic dipoles induced by 
the magnetic templating on the vortex dynamics, it is 
still uncertain whether the ultimate source of disorder in 
the as-grown state is the random distribution of magnetic 
poles in the buckle state or the presence of triangles in 
the magnetic vortex state along one line of connecting 
triangles. In addition, we showed that the high mobility 
of the vortex lattice is a consequence of the generation of 
vortex-antivortex pairs by the magnetic template. This 
picture was confirmed experimentally by high frequency 
transport measurements and low frequency ratchet ex- 
periments. Due to the self organized commensurability 
of these channels a pronounced mode locking effect is ob- 
served. These features were used to probe and tune the 
discrete motion of the composite lattice of vortices and 
antivortices. This combination of techniques represent a 
powerful tool for understanding the motion of vortices at 
microscopic scale. 
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